Violation of Geometrical Scaling in pp Collisions at NA61/SHINE 
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We analyze geometrical scaling (GS) of negative pion multiplicity px distributions at 
NA61/SHINE energies. We show that even though NA61/SHINE energies are low, one may ex- 
pect to find GS in the particle spectra. We argue that qualitative behavior of ratios of multiplicities 
at different energies is in agreement with a simple picture of GS which is violated for px smaller than 
some nonperturbative scale A and when larger Bjorken x of one of the scattering patrons crosses 
i max above which gluonic cloud becomes dilute and quark degrees of freedom become important. 



Scaling properties of physical observables have been 
for a long time one of the most instructive ways to get 
information about the underlying physics. In some cases 
this information could be fully quantified; here the most 
prominent example is Bjorken scaling of Deep Inelastic 
Scattering (DIS) structure functions [T|. In some other 
cases, like the KNO scaling [5] for example, the underly- 
ing physics is still poorly understood. Some theoretically 
conjectured scaling laws, like Feynman scaling in high 
energy hadronic collisions [3], have not been experimen- 
tally confirmed. In any case both the emergence of a 
scaling law, and - sometimes more importantly - its vio- 
lation, test our understanding of microscopic phenomena 
that are of importance in given experimental conditions. 
Again Bjorken scaling is here a benchmark, since both its 
emergence and violation are fully understood in terms of 
perturbative quantum chromodynamics (QCD) [3]. 

In this paper we will focus on another scaling law, 
called geometrical scaling (GS), which has been intro- 
duced in the context of DIS [5] . It has been shown that 
GS scaling is present also in high energy pp collisions 
O [7] . An onset of GS in heavy ion collisions at RHIC 
energies has been reported in Refs. [8]. 

The proton seen at low Bjorken x < ir max is domi- 
nated by a gluon cloud characterized by a typical scale 
Q s (x) which is called saturation scale [9]. Phenomenon 
of saturation (for a review see jTUl E] ) appears due to 
the nonlinearities of parton evolution at small x given by 
so called JIMWLK hierarchy equations [T2] which in the 
large N c limit reduce to the Balitsky-Kovchegov equa- 
tion [13]. These equations admit traveling wave solu- 
tions which explicitly exhibit GS |14j . An effective the- 
ory describing small x regime is Color Glass Condensate 
[T51 [T5] . For the present study the details of saturation 
are not of primary interest, it is the very existence of the 
saturation scale which plays the crucial role. 

Recently it has been shown [17] that GS in DIS works 
very well up to relatively large £ max ~ 0.1 (see also 
[Tg]). This should not come as a total surprise since it is 
known that GS scaling extends well above the saturation 
scale both in the DGLAP [H] and BFKL gU] evolution 
schemes not only for the boundary conditions that ex- 
hibit GS |21j to start with. For x < x max in the first ap- 
proximation the virtual photon cross-section <r 7 * p in DIS 



is a function of a scaling variable r = —q 2 /Q 2 (x) only, 
where q 2 is photon virtuality [S] [5] . Similarly in hadronic 
collisions charged particle multiplicity in central rapidity 
exhibits GS [6] 
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where F is a universal dimensionless function of the scal- 
ing variable 



P \/Ql{x). 
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In this paper we shall address a question of applicabil- 
ity of Eq. ([!]) beyond central rapidity region. This is 
an important and interesting question for the following 
reasons. Firstly, Bjorken x's of scattering gluons that 
produce a particle of a given transverse momentum px 
and rapidity y read 



x h2 = e ±y p T /W 



(3) 



and can be quite different for large rapidities (we shall 
assume positive y, hence x 2 < Xi). Here W — y/s de- 
notes pp scattering energy. If x 2 and x\ are different 
then two different saturation scales emerge: Q s \ and Q s2 
|22) . Whether this leads to violation of GS depends on 
the form of the saturation scale (for discussion of differ- 
ent forms of scaling variable see e.g. 23 and references 
therein). Secondly, and - as we shall argue below - more 
importantly, when the larger Bjorken x = x\ for some 
y leaves the domain of GS, i.e. x\ > x max , violation of 
GS should appear due to the fact that one of the scat- 
tering gluonic clouds becomes dilute and quark degrees 
of freedom become important. We are going to show in 
a model-independent way that this is what indeed does 
happen. 

Theoretical formula for the multiplicity density in pp 
collisions is known since 1981 [24] 

, = ~2 ( d 2 k T a s ipi(xi, k\)ip 2 {x 2 , {k-p)%). (4) 

dyd z p T p^ J 

Here l/p T corresponds to gg — > g cross-section, ipi j2 are 
unintegrated gluon densities and a s stands for strong cou- 
pling constant. Although in principle Q describes gluon 
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production, one uses parton-hadron duality [55] to argue 
that it can be also applied to particle multiplicities where 
constant C takes care of the hardonization effects. 

At high energies and small x's unintegrated gluon den- 
sities depend on fcp through the ratio k\/Ql{x) [5] - 
which is an underlying mechanism behind geometrical 
scaling. Throughout this paper we shall use the follow- 
ing form of the saturation scale [9] : 

Q S 2 K 2 ) = Q 2 (e ±v PT/(xoV~s))- X (5) 

although other possibilities have been also discussed in 
the literature [23.. Parameters Qq ~ 1 GeV/c and 
xq ~ 10~ 3 depend on a particular model of gluon densi- 
ties, whereas exponent A can be extracted from the data 
in a model-independent way. Indeed, by analyzing mul- 
tiplicity spectra at CMS we have found that for pp colli- 
sions A = 0.27 [5]. Note that for A = saturation scale is 
equal to Qq and scaling variable r = p\jQ\ is essentially 
transverse momentum squared in units of Qq. 

There are many sources of possible violation of GS in 
multiplicity spectra given by Eq. (JlJ . Firstly - as already 
discussed - one of the Bjorken x's ([3| can be outside the 
GS domain. Secondly, gluons of very low fcx will feel non- 
perturbative effects of the order of Aqcd- Next, running 
of the strong coupling a s will evidently produce loga- 
rithmic corrections to GS. Finally, gluon fragmentation 
into physical particles will introduce dependence on their 
masses, which should induce GS violation in spectra of 
identified particles. 

In this paper we shall concentrate on the first two 
effects which can be studied by looking at rapidity de- 
pendence of the particle multiplicities in pp collisions. 
For particles produced with y =^ Bjorken x's Q are 
different and we have in fact two different saturation 
scales [52]. Hence instead of F(r) in Eq.Q we will have 
some other function G(ti,t 2 ). This by itself does not 
mean that GS is necessarily violated since we can always 
rewrite G(r 1 ,r 2 ) = G(t 1)( o) with p — Q s i/Q S 2- For fixed 
rapidity y ratio p is constant for Q s of the form given 
by Eq.([5| and universal dependence of multiplicity dis- 
tribution upon scaling variable t\ is retained. Note that 
F(r) = G(r,l). 

On the other hand by looking at the spectra with in- 
creasing y one can eventually reach x% > x max and GS 
violation should be seen. Luckily we have at our disposal 
very recent pp data from NA61/SHINE experiment at 
CERN [26] where particle spectra at different rapidities 
y = 0.1 -j- 3.5 and 5 different energies have been mea- 
sured for p T = 0.025 -7- 1.375 GeV/c. For the present 
analysis we shall use negative pion spectra, which has 
an advantage that one avoids possible GS violation due 
to particle masses. The only problem with this piece 
of data is that the scattering energies are rather low: 
Wi...5 = 17.28, 12.36, 8.77, 7.75 and 6.28 GeV. Nev- 
ertheless, as we will show below, Bjorken x's involved 
are small enough for GS to be present. This is mainly 
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FIG. 1. Kinematic range of NA61/SHINE experiment in 
(W,pt) plane. Data points are shown as stars, squares, cir- 
cles and triangles. Shaded areas denote regions where GS 
is expected to be violated for y — 0.1. For other rapidities 
(marked at the r.h.s. of the plot) GS violation takes place 
above the respective dashed line. 

due to the before-mentioned fact that in DIS £ max ~ 0.1 
[T71 ITS] ; for the purpose of the present analysis we shall 
take x max = 0.09. Then GS should be seen if x\ < x max , 
i.e. for 

PT < PTmaxtW 7 , y) = x maK W e' y . (6) 

On the other hand transverse momenta of produced glu- 
ons should be larger than some nonperturbative scale A: 

PT > PTmin = A (7) 

for which we take A = 250 MeV [27] • Having fixed kine- 
matical range where GS should be present we can have a 
look at kinematical range of NA61 /SHINE spectra shown 
in Fig. [TJ Shaded area below Pxmin = A is excluded for 
all rapidities. Shaded region in an upper part of the 
plot is excluded for y = 0.1, and dashed lines represent 
PTmax(W, y) as functions of W for different rapidities y. 
Therefore for given rapidity only points below these lines 
should exhibit GS. Experimental points for different ener- 
gies are represented by stars (17.28 GeV), squares (12.36 
GeV), circles (8.77 GeV) and triangles (7.75 and 6.28 
GeV). One can see from Fig. [I] that for y — 0.1 GS re- 
gion extends towards the smallest NA61/SHINE energy. 
This is due to the fact that x max is as large as 0.09. By 
increasing y we see that some points fall outside GS win- 
dow and finally for y > 1.7 no GS should be present in 
NA61/SHINE data. One may note that both x max and A 
together with exponent A could be determined in a self- 
consistent way from the NA61/SHINE data. This is be- 
yond the scope of the present note and will be presented 
elsewhere. Therefore our analysis has to be considered 
mostly as a qualitative one. 
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FIG. 2. Ratios as functions of yfr for the lowest rapidity y = 0.1: a) for A = when y/r = pt and b) for A = 0.27 which 
corresponds to GS. 



Limits ([6| and ^ define GS window r m ; n < r < r max 
in scaling variable r 

7"max = ^ ^max^ /^?0 (^max/ ^o) ) 

r m i„ = e^A 2 /Q2 (A/(xoH/)) A^ (8) 

We see from ([8]) that for fixed W, r max falls exponentially 
with rapidity y much faster than r m - m is rising with y. 
When r m j n ~ r max window for GS closes, and we get from 
^ that the larger W the larger y when this happens. For 
fixed y GS window is growing with W. 

Now we come to the formulation of a model- 
independent criterion that allows to asses whether GS 
is present or not. To this end we follow Refs. [BJ [T7] and 
define ratios of multiplicities 

Rik(r;y) = ^-( Wi ,r; y )/-^-(W k ,T;y) (9) 

which, according to , should be equal to unity if GS is 
present. Since GS window is the largest for W\ = 17.28 
GeV, we form 4 ratios Rik with k = 2, . . . 5 where index 
k refers to NA61 /SHINE energies mentioned above. In 
Fig. [2] we plot ratios Rik as functions of y/r for the lowest 
rapidity y = 0.1 firstly for A = 0, i.e. essentially as 
functions of pr, and secondly for A = 0.27. We see that 
for A = all ratios Rik are larger than 1 and rise with r 
up to Rik ~ 1.5-7- 6. In contrast, for A = 0.27 for all ratios 
Rik there exists a domain in y/r where, within the errors, 
they are equal to 1. We see that r m i n for W2 = 12.36 GeV 
is approximately 0.6 while for W$ — 6.28 GeV it is 0.8 
or so. On the other hand for W5 GS windows closes at 
\Jt ~ 1.5, whereas for W2 it extends further than 2.1. 
This is in qualitative agreement with limits ([8]), and as 
we will see in the following, this trend continues once we 
increase y. 

Looking at Fig. [T] we see that there is no qualitative 
difference between rapidities 0.1 and 0.3, and this is what 



indeed can be seen in Fig. |3]a. However, already for 
y = 0.7 the lowest energy data are no longer in the scaling 
region, and we see this clearly in Fig. [3]a. Further on, 
for rapidity y — 0.9 the second lowest energy data leave 
GS window, and this is again nicely substantiated by the 
behavior of R14 shown in Fig. [3]c. Finally for y = 1.3 
none of the tested energies W2 -j- W5 should exhibit GS, 
and this is clearly confirmed by Fig. |3]d. Note that in 
all these cases the reference points at scattering energy 
Wi remain in GS window, therefore scaling violations of 
ratios Rik should be attributed to the points at lower 
energies Wk < Wi- 
lt is somewhat surprising that relatively low energy 
NA61/SHINE data do show GS at all. This is due to 
the fact that GS scaling works up to much higher x's 
than originally anticipated: x max ^0.1. This has been 
confirmed by direct analysis of DIS in Ref. [17] and also 
indirectly by the present study. 

In this note we have explored NA61/SHINE [26 nega- 
tive pion spectra in pp collisions taken at different rapidi- 
ties y. In order to study GS in multiplicity distributions 
we have formed ratios ^ that should be equal to 1 for 
scaling variable r within GS window We have ar- 
gued that behavior of ratios Rik with increasing y is in 
qualitative agreement with the following picture of par- 
ticle production in hadronic collisions. For low pi's GS 
is violated due the existence of the nonperturbative en- 
ergy scale A. For high transverse momenta, larger of two 
Bjorken x's, namely X\, crosses the maximal x for which 
GS is present, Xy > x max , and GS is again violated. For 
rapidities larger than 1.3 the window of GS closes for all 
energies Wk but W\, and no GS is seen in ratios Rik- 
This behavior of GS for different energies with varying y 
is strikingly in qualitative agreement with the kinemat- 
ical constraints shown in FigJT] and more quantitatively 
with the range of GS window given by Eq.(|8]). 

Our aim here was to point out that model-independent 
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FIG. 3. Ratios R\k as functions of \pr for A = 0.27 and for different rapidities a) y — 0.3, b) y = 0.7, c) y = 0.9 and d) t/ = 1.3. 
With increase of rapidity, gradual closure of GS window can be seen. 



analysis of the multiplicity distributions at different ra- 
pidities can provide an interesting insight into the pro- 
duction mechanism of low and medium px pions from the 
point of view of gluon saturation and geometrical scaling. 
One obvious reservation is that for higher rapidities one 
can force ratios R\ k to approach unity at the expense of 
adjusting A. However, very soon A gets bigger than 1 and 
for rapidities of the order of 2 it has to be bigger than 
3. The quality of such "geometrical scaling" is very poor, 
and - more importantly - there is no physical picture 
supporting such big values of A. This is one of the rea- 
sons for which one should perform a more quantitative 
analysis of GS in pp collisions, which will be presented 
elsewhere. 

Let us finish by a remark that experimental data at 
high LHC energies would be best suited for an analysis 
outlined in this letter. For the LHC energies GS closes 
for very high rapidities up to y ~ 8. Since one needs 
data at different rapididities, perhaps the best detector 
for this purpose would be LHC-b. 
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